Artificial electromagnetism for neutral atoms: Escher staircase and Laughlin liquids 
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We show how lasers may create fields which couple to neutral atoms in the same way that the 
electromagnetic fields couple to charged particles. These fields are needed for using neutral atoms as 
an analog quantum computer for simulating the properties of many-body systems of charged particles. 
They allow for seemingly paradoxical geometries, such as a ring where atoms continuously reduce 
their potential energy while moving in a closed path. We propose neutral atom experiments which 
probe quantum Hall effects and the interplay between magnetic fields and periodic potentials. 
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Recently, many researchers have expressed interest in 
using ultracold alkali atoms as analog quantum computers 
to simulate properties of sohd state systems . For ex- 
ample, the leading model of high temperature supercon- 
ductivity, the Hubbard model, can be studied by placing 
alkali atoms in an optical lattice - a periodic potential 
formed by interfering several laser beams. Experimental 
realizations of the Hubbard model could show whether it 
captures the phenomena of high temperature supercon- 
ductivity. Similarly, cold gases provide an ideal setting 
for studying models of quantum Hall effects and exotic 
phase transitions 0- 

A major impediment to studying some of these mod- 
els, such as those describing quantum Hall effects, is the 
lack of fields which couple to the neutral atoms in the 
same way that the electric and magnetic fields couple to 
charged particles. Here, we show how to create these ar- 
tificial electromagnetic fields. Since these fields are only 
analogies of the real electric and magnetic field they do 
not obey Maxwell's equations. One can therefore cre- 
ate unphysical and counterintuitive field configurations 
which lead to a set of as-yet unstudied behavior. Among 
our examples of these seemingly impossible field config- 
urations, we describe an 'Escher staircase' setup where 
atoms can move around a closed path, continually reduc- 
ing their potential energy. 

The literature already contains several, somewhat lim- 
ited, implementation of electrical and magnetic fields 
for neutral atoms. Experimentalists routinely use the 
Earth's gravitational field as an analog of a uniform elec- 
tric field j^]. They also study systems in non-incrtial 
frames: uniform acceleration is equivalent to a constant 
electric field [3, while circular motion corresponds to a 
uniform magnetic field |^. Recently, Jaksch and ZoUer 
described a method where an effective magnetic field 
can be applied to two-state atoms in an appropriately de- 
signed optical lattice in the presence of an external 'elec- 
tric field'. Our approach is an elaboration of Jaksch's, 
where the two-state atoms are replaced by three-state 
atoms. This allows us to overcome the major deficiency 
of Jaksch and Zoller's scheme: we do not need an ex- 
ternal electric field to generate the magnetic field. This 
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FIG. 1: (Color Online) Lattices of three types of sites: (a) 
1-D chain, (b) square, (c) ring, and (d) triangular. 



improvement comes at the cost of more complicated laser 
configurations. 

As in these prior analogs of electromagnetism, our ar- 
tificial fields contain no dynamical degrees of freedom. 
Therefore they neither give rise to analogs of Coulomb 
interactions between the neutral atoms, nor do they sup- 
port analogs of light. 

Subsequent to the preparation of this manuscript, an- 
other scheme for generating analogs of electromagnetic 
fields was suggested by Sorensen, Demler, and Lukin 
0. That work uses time dependent hopping matrix ele- 
ments along with a large oscillating quadrupolar poten- 
tial. Compared to our approach, Sorensen et al. use a 
much simpler laser configuration, however there are non- 
trivial technical issues involved with the stability of the 
oscillating potential. 

Basic Setup: Our approach relies upon creating an 
optical lattice with three distinct sets of minima. Each 
of these minima trap a different internal state of the neu- 
tral atoms. The internal states will be labeled 'A', 'B', 
and 'C, and the minima will be labeled by their location 
and by the state that is trapped at that location. For ex- 
ample, figure 1(a) shows a one-dimensional array labeled 
as • • -Ai-B2-C3-A4-B5-C6-A7-B8-- • •. Importantly, this 
setting breaks parity symmetry. 

Looking at this one dimensional chain, an atom in state 
A, sitting in site A4, is immobile. The atom cannot hop 
to site B5 or C3, because it would need some mechanism 
for changing its internal state. The probability of tun- 
neling by three sites to Ai or A7 is astronomically small. 

We turn on hopping between site A4 and B5 by in- 
troducing a laser with the following properties: (i) the 
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laser frequency ujab is close to the energy differenced be- 
tween the internal states A and B (ie. luab ~ Ea — Eb); 
(ii) the laser polarization is chosen so that the transition 
from internal state A to B is allowed; (iii) the laser can- 
not induce transitions from states A to C or from B to 
C, either because the transition is forbidden, or because 
the detuning is too great. One does not have to use a 
single laser to drive this transition but can instead use a 
Raman transition, which involves multiple lasers and the 
virtual occupation of one or more intermediate state. 

In the presence of this laser field, the atom can explore 
a two state Hilbert space. In the rotating wave approxi- 
mation, the time dependent Schroedinger equation is 

_( Ea -f7ABe-*("*^*+"*) A 

^^^'^ [-nABe'^^^^'+^^ Eb )' 

The quantum mechanical amplitude for the particle being 
in state A (B) on site A4 (B5) is -0^14 ("035)- The energy 
of the internal states A/B are Ea/b- The Rabi frequency 
f^AB is proportional to the product of the laser amplitude 
and the overlap between the states trapped in A4 and B5 . 
We take ^Iab to be real, and introduce a phase 4>, which 
is related to the phase of the coupling laser. In particular, 
if we translated the entire lattice by some distance r, the 
phase would change by ^ -I- q • r, where q is the 
wave- vector of the coupling laser . 

This, and future Hamiltonians arc more compactly 
written in a second quantized notation, 

H = EAi'i4i^A4 + EBlplr,TpB5 (2) 

where, for example, creation and annihilation operators 
^]^^ and '0A4 add and remove a particle from site A4 in 
internal state A. In the non-interacting system, the oper- 
ators obey the same equations of motion as the wave- 
function 0; in (^). At the single-particle level it does not 
matter whether we use bosonic or fermionic commuta- 
tion relations. Where no confusion will result, we may 
neglect the letter A which denotes the internal state. 

We apply time-dependent canonical transformations 
of the form 0^ e^^^^^'ipj, -01 e'^^^^^ijj^y As is 
readily verified from the equations of motion un- 
der this transformation the Hamiltonian becomes H 
H — f'{t)ipjtl;j. In particular we can construct a time 
independent Hamiltonian by transforming into the 'ro- 
tating frame,' 

0'A4 - e-*(^**-«^A4 (3) 
^B5 ^ e-^(^Bt+A.B)^gg (4) 

H = ~T(^1^5+V3^V'4) + A7At^5, (5) 

where r = ^ab and A = lo ab ~ (-^a ~ Eb)- 



Introducing two more lasers, coupling states B-C, and 
C-A with appropriately chosen intensities and detunings, 
this same procedure yields the Hamiltonian 

H = (j-A(0]0,) - t(0]v3,+i -t- V-l+i^,)) , (6) 

corresponding to a one-dimensional chain of sites in a 
uniform electric field. As is shown below, this same ap- 
proach can produce electric fields in higher dimensions. 
In this case, momentum transfer from the lasers will gen- 
erate an effective magnetic field. 

Higher Dimensions: In more complicated geome- 
tries there may not be a Canonical transformation which 
leads to a time independent Hamiltonian. However, the 
time dependence takes a simple form if one transforms 

where j labels the site located at rj, and /ij = A,B,C 
gives the internal state which is trapped at that site. The 
Hamiltonian then becomes 

H = -Y.^,.H (e^'^-^^-^'^e-'^^--^ Vi^.V-^,, + H.C.) . 

(8) 

The sum includes all nearest neighbor sites (ij). The 
internal state trapped at site i is /i^. The bond position 
is Rij = (ri -I- rj)/2. The hopping is r^^ = fi^,^ for 
and r^^ — tq. The parameter tq is given by the overlap of 
the wavefunctions on neighboring sites. The wave-vector 
of the laser couphng state /i to is q^^ (so q^^ — 0). 
The detuning is A^j^ = uj^i, — {E^ — E^) when ji ^ v, 
and A^^ = 0. The letters H.C. denote the Hermitian 
conjugate of the previous term. 

If all of the laser intensities are adjusted so that t^i, — 
To for all /i, t^, then equation is formally the equation 
of motion of a particle with charge e in a vector potential 
defined on the bonds by 

^A(Rij) • rij = q^,^^. • Rij - A^^^^t, (9) 

where rij = ri — rj . 

Using this mapping to a vector potential, we can con- 
struct many interesting field configurations. For exam- 
ple, consider a lattice with the striped geometry shown 
in figure 1(b), where as one moves in the x direction, one 
encounters alternating rows of sites A, B, and C. With 
this geometry, only the x-component of the vector po- 
tential. Ax, will be non-zero. In the simplest case, where 
each of the three coupling lasers have the same wave- 
vector q and detuning A, the vector potential is A(r) = 
x(c/ec?)(q • r — At), where d is the lattice spacing. This 
corresponds to a uniform electric field E = — xAc/ed and 
a uniform magnetic field B = x x q(c/ed). By changing 
the relative angle between q and the x axis, one can con- 
trol the strength of the magnetic field. Since the recoil 
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momentum q can be made comparable to the inverse lat- 
tice spacing, one should be able to construct extremely 
large fields where flux through a unit cell of the lattice 
exceeds the magnetic flux quantum $0 = 27rc/e. 

If q is aligned with the hopping direction, then the 
effective magnetic field vanishes, resulting in an electric 
field without a magnetic field. 

Applications: Earlier we introduced some interesting 
problems which could be addressed by applying effective 
electric and magnetic fields to a system of particles on a 
lattice. Here we discuss a further possibilities. 

At moderate values of the "magnetic field" exper- 
iments could explore how the periodic potential af- 
fects vortex structures in a Bose condensate |0. One 
could also study vortex physics near "pairing transitions" 
where the structure of vortices change 0| 

At much larger fields (<I> ~ $0) Jaksch and ZoUer re- 
cently discussed the exciting idea of using neutral atoms 
to study the fractal energy spectrum that Hofstadter 
predicted for noninteracting charged particles on a lat- 
tice in a magnetic field. The spectral gaps would be 
observable as plateaus in the density of noninteracting 
harmonically trapped fermions. It would be even more 
exciting to explore an interacting system in this same 
regime, and study fractional quantum Hall physics, and 
the interplay between quantum Hall effects, Mott insu- 
lating physics, and this fractal single-particle spectrum 
[l3j . The simplest such experiment would use the geom- 
etry in figure 1(b), and allow the system to equilibrate 
with A = 0. All single particle observables are measur- 
able through imaging, while photoassociation provides 
access to the short range pair correlation function^]. 
Some transport measurements are possible by detuning 
the lasers so that A 0. 

Several authors have shown that for filling fractions 
1/2 < V < 6, bosons with short range interactions in 
a strong magnetic field will form non-trivial many-body 
states Fermions are more tricky, as s-wave interac- 
tions (which dominate at low temperatures) cannot lead 
to fractional quantum Hall effects in fermions. However, 
resonantly enhanced p-wave interactions can lead to such 
correlated states [l^ . 

Previous proposals for creating analogs of quantum 
Hall states in cold atoms relied upon rotation to provide 
the effective vector potential. Such schemes are made 
difficult by the need to carefully balance the centripetal 
force which maintains rotation and the harmonic trap- 
ping potential. The window of rotation speeds for find- 
ing strongly-correlated physics falls off with the inverse 
of the number of particles. The present approach does 
not require this delicate balancing of forces, and there- 
fore allows one to study these effects in a macroscopic 
system. 

Not only are magnetic fields of interest, but so are 
large electric fields. For example Sachdev et al. 0| have 
discussed the intricate Mott-Insulator states which are 



found when the 'voltage difference' between neighboring 
wells is comparable to the on-site repulsion. The method 
presented here is a powerful tool for studying such states. 

Unphysical Fields: We once again emphasize that 
although A couples to the neutral atoms as if it were 
a vector potential, it does not obey Maxwell's equa- 
tions. Consequently, one can engineer seemingly para- 
doxical geometries. Consider, for instance, the ring of 
sites illustrated in figure Ic, with all detunings set equal. 
According to equation Q, there is a uniform 'electric' 
field pointing along the chain. Thus a particle can move 
around the ring, continuously moving to a lower poten- 
tial energy, returning to the starting point, but (by con- 
servation of energy) having gained a great deal of ki- 
netic energy. One can repeat the process ad infinitum; 
the maximum velocity is limited only by Umklapp pro- 
cesses. That is, when the particles deBroglie wavevector 
coincides with the intersite distance, the matter-wave is 
Bragg refiected off of the lattice, and reverses direction. 
If the chain was not bent in a circle, this reflection would 
lead to the familiar Bloch oscillations. No conservation 
laws are violated by this continuous acceleration, as the 
lasers provide a source of energy and momentum. 

This bizarre situation where a particle can reduce its 
potential energy by moving in a closed path is rem- 
iniscent of the optical illusion in MC Escher's print 
"Ascending and Descending," where a staircase forms 
a continuously descending closed loop. The quantum 
mechanical properties of a particle in such a chain of 
N sites are ascertained by noting that the Hamilto- 
nian, H = -T^f^i(e^**^]V,+i + e-^^^Vj+i^j), with 
TpN+i = ^1, is translationally invariant, and therefore 
extraordinarily simple in momentum space. In terms 
of operators — e^'^^^^^^'^ tjjj / \/N , the Hamilto- 

nian is diagonal, H = '^f. Ek(t)a\.ak. The eigenvalues 
Ek{t) — — 2r cos(27rfc/A + St) are time dependent, re- 
flecting the non-equilibrium nature of the system. The 
motion of a wave-packet is determined by the instanta- 
neous phase velocity 



= 2Tdsm{2TTk/N + St), 



(10) 



which oscillates as a function of time. The factor of 
dN/2'K, where d is the intersite spacing, converts the 
velocity into physical units. This oscillation is exactly 
the Bragg diffraction previously mentioned. During one 
period of oscillation, the particle moves around the ring 
approximately 2t / {N5) times. 

A more complicated geometry with similar paradoxical 
properties is illustrated in figure Id. In this structure, a 
triangular lattice is formed from three interpenetrating 
sublattices with wells of type A, B, and C. Here, a con- 
stant detuning yields a very intricate 'unphysical' elec- 
tric field configuration: arrows depict directions in which 
hopping reduces the potential energy. Upon traversing 
alternate plaquettes, a particle can continuously increase, 
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or decrease its potential energy. To understand the be- 
havior of a particle in this lattice, one once again relies 
upon translational invariance, and introduces operators 
Ofe = ^rS"*'^ '^, where k lies in the first Brillouin zone 
(BZ) of the triangular lattice, and the sum is over all 
lattice sites. The Hamiltonian is then 

H = -tJ2 (E-=iV'JV'.+,,) +H.C.] (11) 

r 

= -2ry^_^^4afeE-=iCos(k.rj+5t). (12) 

The lattice generators {ri, r2, ra} connect nearest neigh- 
bor sites, and are illustrated by arrows in figure Id. 
Only two of these generators are linearly independent 
(ri -I- r2 -I- rs = 0). The area of the first Brillouin zone is 
ri = Stt'^/ VSd^, where d is the lattice spacing. Again, the 
group velocity of a wave packet is simply the gradient of 
the energy Ek = — 2t^^- cos(k • rj -|- 6t). Of particular 
note is the fact that at the zone center (fc — 0) the group 
velocity is alway zero. Thus a stationary packet remains 
stationary. This result is not surprising, since there is 
nothing in the geometry which picks out a direction in 
which the packet could start to move. 

More surprising is the fact that the effective mass, re- 
lated to the curvature of Ek is oscillatory at fc = 0, spend- 
ing equal amounts of time positive and negative. When 
the effective mass is negative, quantum diffusion acts op- 
posite to its normal behavior, and wave packets become 
sharper. Thus localization occurs: the wave packet's 
size oscillates periodically, rather than continually grow- 
ing. Similarly, if the packet has a small momentum with 
I A; I <C 27r/d, then the particle does not simply propagate 
ballistically, but its velocity oscillates sinusoidally about 
?; = 0, and the particle is trapped near its initial location. 

Physical Realization: There are many ways to en- 
gineer the three-state lattices described above. The diffi- 
cult task is to produce the confinement and Raman cou- 
plings with a small number of lasers in a geometry which 
can be easily implemented. A detailed analysis of the var- 
ious configurations goes beyond the scope of this paper, 
and a more comprehensive article is in preparation. 

A key idea is that if the internal states are related by 
symmetries (ex. a spin-1 multiplet), then the various 
traps can be created by the same lasers, and the {A-B) 
and {B-C) Raman transitions can use the same drive. 
Driving transitions with microwave or RF fields, rather 
than lasers, will reduce the need for optical access [l^ . 

An alternative approach is to note that one can create 
analogs of electromagnetic fields even if the sites A, B, 
and C, trap atoms in the same state. One can instead 
rely on a superlattice structure, where the energies of the 
three sites differ by large amounts [l7| . Hopping is only 
possible if a Raman laser supplies the missing energy; 
detuning and recoil give the same effects as in the case 
with different internal states. 
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